Modifying Fragility and Collective Motion in Polymer Melts with Nanoparticles 
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We investigate the impact of nanoparticles (NP) on the fragility and cooperative string-like motion 
in a model glass-forming polymer melt by molecular dynamics simulation. The NP cause significant 
changes to both the fragility and the average length of string-like motion, where the effect depends 
on the NP-polymer interaction and the NP concentration. We interpret these changes via the 
Adam-Gibbs (AG) theory, assuming the strings can be identified with the "cooperatively rearranging 
regions" of AG. Our findings indicate fragility is primarily a measure of the temperature dependence 
of the cooperativity of molecular motion. 



The addition of a small concentration of nanoparti- 
cles (NP) to glass-forming (GF) polymer materials can 
lead to large property changes that are difficult to com- 
prehend by extension of the effects of macroscopic filler 
additives. Depending on system details, changes may be 
rationalized by the large surface-to-volume ratio of NP, 
chain bridging [TH3], or NP self-assembly into extended 
structures [4j [5] . Changes in the glass transition temper- 
ature Tg have been emphasized, and both experimental 
and theoretical studies indicate that attractive or repul- 
sive (non-attractive) polymer-NP interactions tend to in- 
crease or decrease Tg, respectively. Correspondingly, the 
interfacial polymer layer around the NP shows a slowing 
down (increased Tg) or an acceleration of dynamics (de- 
creased Tg), providing a molecular interpretation [BHS]. 

Unfortunately, Tg changes provide only a limited un- 
derstanding of how NP affect the properties of GF poly- 
mer melts. We also expect that the T dependence of 
dynamical properties approaching Tg, or the "fragility" 
of glass formation [TD], will be altered. Fragility changes 
have been argued for on theoretical grounds |TT], based 
on the finding that changes in the molecular packing in 
the glass state [T < Tg) should also alter the fragility 
of glass formation. The NP we study should be particu- 
larly effective at modifying molecular scale packing and 
motion, since their size is roughly commensurate with 
the heterogeneity scale of fluids near Tg [T^ . 

In this letter, we address how polymer-NP interactions 
and NP concentration affect the fragility of glass forma- 
tion, and how this effect relates to cooperative motion. 
We demonstrate that the changes in the relaxation time 
r can be related to changes in the average size L of the 
string-like cooperative motion of monomers. The Adam- 
Gibbs (AG) theory [l^], predicts a specific relationship 
between between size of hypothetical "cooperatively re- 
arranging regions" (CRR) and structural relaxation time. 
If L corresponds the size of the CRR, we find that the 
AG relation holds for all concentrations </> and interaction 
types considered. Additionally, the AG theory predicts 
that fragility is sensitive to the T dependence of the size 



of CRR, and we confirm this relationship. 

Our findings are based on equilibrium molecular dy- 
namics simulations of a nanoparticle surrounded by a 
dense polymer melt, as well as simulations of a pure melt 
for comparison purposes. We utilize periodic boundary 
conditions so that our results represent an ideal, uni- 
form dispersion of NP. The polymers are modeled by a 
well-studied bead-spring model [H], but with the cut- 
off distance between pairs extended to include attractive 
Lennard- Jones (LJ) interactions. All monomer pairs in- 
teract via a LJ potential, and bonded monomers along a 
chain are connected via a FENE anharmonic spring po- 
tential. The NP consists of 356 Lennard- Jones particles 
bonded to form an icosahedral NP; the facet size of the 
NP roughly equals the equilibrium end-to-end distance 
for a chain of 20 monomers. Details of the simulation 
protocol and our model potentials can be found in the 
supplemental material and in ref. [5] . 

We simulate systems with 100, 200, or 400 chains of 
M = 20 monomers each (for totals of = 2000, 4000, 
and 8000 monomers) to address the effect of varying the 
NP volume fraction. Under constant pressure conditions, 
the addition of nanoparticles can give rise to a change in 
the overall melt density. A slight change in density can 
cause a significant change in the dynamic properties rel- 
ative to the pure melt. In order to probe only changes 
caused by the interactions between the NP and the poly- 
mer melt, we have matched the density of monomers far 
from the NP with that of the pure polymer melt [5]. 

To quantify changes in the nanocomposite dynamics, 
we evaluate the effect of and the polymer-NP interac- 
tions on r, measured from the relaxation of the coherent 
intermediate scattering function (see supplementary in- 
formation). The effects of interactions on r and Tg for 
some were presented in ref. [5]; here we provide addi- 
tional simulation data and focus our analysis on fragility 
and cooperative motion. As expected. Fig. [T] shows that 
attractive polymer-NP interactions slow the relaxation 
(t becomes larger), while non-attractive polymer-NP in- 
teractions give rise to an increased rate of relaxation (r 
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FIG. 1: Structural relaxation time r as a function of T for 
each (j) normalized by Tpme for the pure melt. The inset shows 
the corresponding Tg as a function of concentration 0. At- 
tractive interactions increase r and Tg, while non-attractive 
interactions decrease r and Tg. In both cases, the effect is 
more pronounced with increasing concentration. The concen- 
trations are (p = 0.0426 (red Q s-nd orange A), 0.0817 (green 
□ and brown <) and 0.151 (blue O and violet v)- The pure 
melt is indicated in black. The fits of the VFT relation to the 
data deviate by at most 0.5 %. 
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becomes smaller). The effect of (p is more clearly seen 
by rescaling r by the value Tpmc in the pure melt, which 
shows that r can be altered by a factor of more than an 
order of magnitude on cooling. The effect of the NP is 
more pronounced at low T. 

We next examine how these changes in r affect Tg and 
fragility. For reference, the inset of Fig. [T] confirms that 
Tg increases when there is attraction, and decreases with 
non-attractive interactions. To estimate Tg, we fit the 
data using the Vogel-Fulcher-Tammann (VFT) expres- 
sion TO] T cx exp {D/[T/Tq — 1]). To is an extrapolated 
divergence T of t, while D provides one measure of the 
fragility. We use the VFT fit to estimate Tg based on 
the condition that T{Tg) = 100 s (the canonical defini- 
tion of the laboratory glass transition [10]), assuming the 
one time unit in standard LJ reduced units corresponds 
to 1 ps (reduced units defined in supplementary informa- 
tion). 

Since there is no single agreed upon measure of 
fragility, we consider several different measures to ensure 
consistency. First, as indicated above, the parameter D 
from a VFT fit to r is widely utilized; specifically, a larger 
value of D indicates a stronger (less fragile) GF fluid so 
that increases with increasing fragility. The most 
common definition of fragility is based on the T depen- 



FIG. 2: (a) Fragility dependence on concentration rela- 
tive to the pure melt. We consider five different measures 
of fragility, which are discussed in the text. All measures of 
fragility show the same qualitative trend: namely, the system 
with attractive polymer-NP interactions becomes more frag- 
ile, while the system with non-attractive interactions becomes 
less fragile (stronger). The last measure, related to the string 
length L, is discussed later in the text, (b) Demonstration of 
the proportionality between Tg and fragility m, as well as the 
high-T activation energy Eoo ■ The different colors of symbols 
indicate the nature of the polymer-NP interactions. 



dence of r near Tg, namely 
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For strong GF systems, the rate of change of r with re- 
spect to T is smaller than that of fragile systems; hence 
m is larger for more fragile GF fiuids. We estimate m 
using our VFT fit. Fragility can also be estimated by 
the ratios To/Tg or Tg/Tc- We estimate Tc using the 
power-law form r ~ (T/Tf. — in an appropriate T 
range (see supplementary information for fitting details) . 
To/Tg and Tg/Tc are larger in more fragile systems [TT] . 

We summarize the results for the various fragility met- 
rics in Fig. |2ja) , where we find that - for all definitions 
- attractive polymer-NP interactions lead to more frag- 
ile glass formation as a function of 4>i conversely, non- 
attractive polymer-NP interactions lead to stronger glass 
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formation. Tliese changes in fragility mirror the changes 
in Tg (fig. [l] inset). In particular, fig. [2]jb) shows that 
m Tg, consistent with experimental trends in pure 
polymeric glass formers [TB] and analytic calculations 
based on the entropy theory of glass formation [T^ . Fig- 
ure [sjb) also shows that the high-T activation energy 
Eoo is roughly proportional to Tg. We discuss the impli- 
cations of these scaling relationships below. 

Our findings for fragility changes are consistent with 
experimental studies of polymer-NP systems. Bansal et 
al. [5] found that dispersions of NP having repulsive inter- 
actions caused Tg to decrease, accompanied by an appre- 
ciable broadening of the glass transition region, indicative 
of increased strength (decreased fragility) of glass forma- 
tion. For fuUerenes dispersed in polystyrene, Cabral and 
co-workers [TH] reported behavior expected for attractive 
polymer-NP interactions, namely an increase in Tg, ac- 
companied by an increased fragility. For small 0, negligi- 
ble changes in the fragility have been reported [12], also 
consistent with our small results. Our results are likely 
not applicable when the NP-polymer interactions are so 
strong that non-equilibrium effects (leading to "bound" 
polymer) [2^ or phase separation may dominate. 

We next examine how the NP interactions impact the 
heterogeneity of molecular motions and how this relates 
to the observed changes in Tg and fragility. Both small- 
molecule and polymeric liquids exhibit pronounced spa- 
tial correlations in mobility, commonly referred to as "dy- 
namical heterogeneity" [T^]. In particular, the most mo- 
bile atoms or molecules tend the cluster on a time scale 
after the "breaking of the cage" , but before the primary 
relaxation r. These clusters of mobile molecules can be 
further dissected into "strings" involving particles mov- 
ing roughly co- linearly |21| ; these structures appear to be 
the most basic units of cooperative relaxation. The char- 
acteristic size of both the mobile-particle clusters and 
strings grow as a fluid is cooled toward Tg. Notably, 
the string-like collective motion is not strongly correlated 
with chain connectivity [Hj , so it should not be confused 
with reptative motion. 

We evaluate the string size L{T) following the proce- 
dures developed in ref. [2Tj . which are slightly modified 
for the specific case of the bead-spring polymer we study, 
as discussed in ref. [22] (see supplementary material for 
the technical details). Figure [s] shows L for all T and 0, 
where we compare systems by normalizing by the value 
of L of the pure melt. Given that L grows as r grows 
on cooling, we expect that the attractive NP interactions 
(which increase r relative to the pure melt) should cause 
L to increase relative to the pure melt, and vice- versa for 
non-attractive NP interactions. Indeed, the variation of 
L is consistent with these expectations, and we conclude 
that the the attractive NP interactions cause an increase 
in the degree of correlated molecular motion for fixed T, 
while non-attractive NP interactions cause a decrease in 
L. Note that the T dependence of L is qualitatively sim- 
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FIG. 3: T dependence of L for different (j) normalized by L 
for the pure melt. Note the parallelism to fig. [T] The inset 
also shows a snapshot of an example string of length L = 9. 



ilar in all cases. 

While the changes in L reflect the changes in r for 
any given T, how do the T dependence of these changes 
compare? In other words, can L be used to predict the 
fragility changes in fig.[2]ja)? To answer this question, we 
are guided by AG theory, which proposes that relaxation 
in GF liquids is dominated by cooperatively rearranging 
regions (CRR). Specifically, AG argue that r is related 
to the average number of particles in the CRR z by a 
generalized Arrhenius relation |23) , 



T = Too ex.p{zEoc/T). 



(2) 



For high T, cooperativity is minimal, so z « 1, and E'oo 
(a constant) can be identified with the high-T activation 
energy of the fluid. Accordingly, the T-dependent activa- 
tion energy E{T) — z{T) Eoo should govern the fragility. 

The strings are a natural candidate to describe the ab- 
stract CRR of AG. Fig. |4] shows that eq. ^ with z (x L 
provides an excellent prediction for r, consistent with 
identifying L with the CRR of AG [22. Therefore, the 
non- Arrhenius behavior of r can be viewed as a conse- 
quence of the increase in L on cooling. Accordingly, L 
must encode the fragility of glass formation. In particu- 
lar, combining eqs. ([T]) and ^ implies a direct relation 
between m and L, 



(Eoo/Tg) [L{Tg)-TgdL/dT\T^ 



(3) 



Since E^o and Tg are proportional (fig. [2jb)), the prefac- 
tor Eoo/Tg is roughly constant for all 4>, and thus does not 
impact the nanocomposite fragility relative to the pure 
melt. Therefore, the fragility changes measured by m on 
adding NP should result from changes to L and T dL/dT 
at Tg. To test this, we extrapolate L to Tg by assum- 
ing consistency between eq. ^ and the VFT expression. 
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FIG. 4: Test of the AG relation (eq. between the size of 
CRR and the structural relaxation time r, assuming z cc L. 
The main figure shows the data scaled by the fit parameters 
so that all data sets collapse to a master curve. Too is defined 
by eq. [2j replacing z hy L. The inset shows that r can be 
collapsed by plotting as a function of T — Tg, a result of the 
proportionality of Tg and m (see fig. [2|b)). 

Figure[2][a) shows that L — TgdL/dT indeed accounts well 
for the observed changes in m. In particular, the contri- 
bution from Tg dL/dT is 3 to 8 times larger than L, which 
ranges from 4.5 to 6.5 at Tg. The range for L is quan- 
titatively consistent with ref. and qualitatively con- 
sistent with experimental evidence for weak sensitivity of 
fragility to the scale of cooperativity |25j[26]. Hence, near 
Tg, fragility is primarily controlled by dL/dT, rather than 
L, consistent with ref. [55] ■ Stated more plainly, fragility 
is primarily a measure of the rate of change dL/dT of the 
extent of cooperative motion. 

The approximate proportionality between m and Tg 
observed in in our system (Fig. [2jb)), and for many high 
molecular mass polymer polymer materials [1^], implies 
an important simplification. Specifically, since m ~ 1/D, 
the product DTg should be nearly constant, and hence r 
should be a nearly universal function of {T ^ Tg). The 
inset to Fig. [3] shows that such a shift collapses all our 
nanocomposite data rather well; this alternate data re- 
duction also indicates the consistency of our Tg extrap- 
olation. However, we are careful to point out that the 
proportionality m cx Tg is not universal |27j . and this 
trend can even be reversed in polyelectrolyte materials 
of interest in battery applications [T71 US] ■ 

In summary, the addition of NP to a polymer melt can 
lead to significant changes in both Tg and fragility, which 
can be related to the cooperative string-like motion. Our 
results support the identification of the strings with the 
abstract CRR of the AG theory, and complement tests 
of the entropy formulation of the AG theory [50] . 
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